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It has recently been suggested that the baryon washout problem of standard elec- 
troweak baryogenesis could be avoided if inflation ends at a low enough energy den- 
sity and a parametric resonance transfers its energy repidly into the standard model 
fields. We present preliminary results of numerical simulations in a SU(2)xU(l) 
gauge-Higgs model in which this process was studied. 

DAMTP-2000-126 
1 Introduction 

In order to explain the observed baryon asymmetry of the universe, a theory 
must satisfy three conditions!!! 

1. Baryon number violation, 

2. C and CP violation and 

3. Deviation from thermal equilibrium. 

In principle, they all seem to fulfilledin the electroweak theory and the stan- 
dard cosmological Big Bang scenario.0 At high temperatures, baryon number 
is violated by non-perturbative sphaleron processes, which change the Chcrn- 
Simons number and consequently, due to a quantum anomaly, also the baryon 
number. The electroweak phase transition makes the system fall out of equi- 
librium in a natural way. The strength of CP violation in the electroweak 
theory is too small, but even that is not a severe problem because the con- 
straints for CP violation arising from beyond the standard model are fairly 
weak. 
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However, the baryon asymmetry generated in the electroweak phase tran- 
sition gets easily washed out. Although sphaleron processes become less fre- 
quent after the phase transition, they don't stop completely. Instead, their 
rate is proportional to exp(— M sp h/T), where M sp h is proportional to the ex- 
pectation value 4> °f the Higgs field, and unless 4> is large enough, the baryon 
asymmetry is washed out. This can only be avoided if the transition is strongly 
enough first order so that the discontinuity of the Higgs field isB A<f> > T. 

In the minimal standard mode], the Higgs mass m# is the only unknown 
parameter, and lattice simulationsu have revealed that, whatever its value, the 
transition is not strong enough. In more complicated models, such as MSSM, 
there are more unknown parameters and this freedom makes it possible to 
satisfy the constraint, but only barely. nn 

In an alternative scenario proposed recently by two groupsJj'EI the baryon 
asymmetry is generated by sphaleron processes during a period of preheating 
after inflation. This requires that inflation ends at an energy scale that is be- 
low the electroweak scale and that a large fraction of the energy of the inflaton 
is transferred rapidly to the standard model fields by a parametric resonance!] 
In the resulting non-equilibrium power spectrum, all the fermionic fields and 
the short-wavelength modes of the bosons are practically in vacuum, but the 
long-wavelength bosonic modes have a high energy density. The sphaleron 
rate depends strongly on the temperature of these long-wavelength modes 
and is therefore very high, and the out-of-equilibrium processes can generate 
a large baryon asymmetry very quickly. Eventually, the system equilibrates 
and the effective temperature decreases by a rate that is much faster than 
the decay rate of baryons. The final temperature Trchoat is determined by the 
initial energy density and provided that it is low enough, T re hcat < 0.6T C , the 
sphaleron rate becomes negligible and the baryon washout is avoided. 

Simulations of the dynamics of preheating in an Abelian gauge-Higgs 
systemtj has confirmed this qualitative picture, but they cannot address the 
issue of baryogenesis directly. In this talk, we discuss the electroweak theory 
with the full gauge group SU(2) xU(l), and argue that reliable simulations 
are possible using reasonable approximations. We present some preliminary 
results and discuss the prospect of determining the generated baryon asym- 
metry using these simulations. 

2 Simulations 

Let us consider a simple model of inflation, in which the expansion of the 
universe is driven by the potential energy of a scalar field, the inflaton, rolling 
slowly down its potential towards its minimum at the origin. Inflation dilutes 
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away all inhomogeneities and thus all the standard model fields are in vacuum 
when inflation ends, and the inflaton has a large homogeneous expectation 
value. It goes on rolling down its potential and starts to oscillate about the 
minimum of its potential. We assume that it is coupled to the Higgs field 
and the two fields start to resonate,Q whereby a large amount of energy is 
rapidly transferred from the inflaton to the long-wavelength modes of the 
Higgs. The details of this process of preheating depend on the properties of 
the inflaton, which are unfortunately unknown. Therefore we simply assume 
that the energy transfer is extremely efficient and results in a state in which the 
energy is concentrated in the Higgs modes with very long wavelengths. From 
the point of view of the microscopic physics that we want to describe, this is 
practically equivalent to the Higgs having a very large value <pQ. Furthermore, 
we assume that the effect of the inflaton to the later dynamics of the system is 
negligible, and therefore we don't include it in our simulations as a dynamical 
field. Thus we can simply consider the time evolution of the standard model 
fields with the special initial conditions in which the Higgs has initial value 
(f>o and all the other fields amin vacuum. 

As previous simulationsatj have shown, the system will quickly reach a 
quasi-equilibrium state in which the long- wavelength modes of the Higgs and 
gauge fields are approximately in equilibrium at a high effective temperature. 
The decays of the bosons transfer energy slowly into the fermions and gluons, 
and the effective temperature of the long- wavelength bosonic modes decreases. 
As the fermions and gluons are initially in vacuum, this process can be de- 
scribed perturbatively by a damping term, whose magnitude T w 2 GeV is 
obtained from the observed lifetime of W and Z bosons. This approximation 
is valid until t ~ T -1 . 

Thus, the only fields that we have to consider are the gauge bosons and 
the Higgs field. For baryogenesis, the relevant degrees of freedom are the 
long-wavelength modes, and they have large occupation numbers. Therefore 
they behave classically, and thus we can study the dynamics of the system 
simply by solving the classical equations of motion 

&l<t> = dm + 2A Q v 2 - ^ 4, - ra o 0, 

dlB, = -djBij + g'lmtf D i( f> - Td B t , 

a$Wi = -[DjtWij] + ig U(Di^ - i(A«A)V " h-c-J - Fd Q W h (1) 
where <fr is the Higgs field, Bi is the U(l) gauge field and Wi is the SU(2) 
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gauge field, and the covariant derivative is 

D i = d i - l -gW i - l -g'B i . (2) 

In addition, both gauge fields must satisfy the corresponding Gauss laws 
d l E i = g'lmir^cf), 

[DuFi] =ig (V - 1^ _ h.c.) , (3) 

where 

TT = d & Ei = -d Bi, Fi = -d Wi. (4) 

Although the classical equations of motion describe the dynamics of the 
long-wavelength modes, they fail to describe the early stages of the thermal- 
ization, when the quantum fluctuations play an important role. Therefore we 
approximate them by adding to the initial field configuration Gaussian fluctu- 
ations with the same two-point correlation function as in the quantum theory 
at tree level. For each real field component Q of mass m and its canonical 
momentum P, this means 

(Q*(t,k)Q(t,k')) = j=L (27r)V 3 >(fc-fc'), (5) 

2Vfc 2 + m 2 

(P*(t,k)P(t,k')) = v/g2 2 +m2 (2^ 3 )(fc-fc0- (6) 

In a sense, this means that the quantum effects arc approximated to lead- 
ing order in perturbation theory Just like real quantum fluctuations, these 
fluctuations generate radiative corrections to the couplings, and they must 
be taken into account, i.e., the parameters must be renormalized. The situ- 
ation is made more complicated by the damping term T, which damps also 
the fluctuations and therefore the mass divergence and consequently the mass 
counterterm also decrease with time. The parameter with the largest ra- 
diative corrections in the Higgs mass, and we have calculated the necessary 
renormalization counterterm at one-loop level in perturbation theory 

2 2 ( 9 o 3 ,,\ 0.226 rt 

Watt « rn% - [6X + - 9 2 + - 9 ' 2 J ^e" r *. (7) 

In a similar way, when we plot {$4>), the quantity is actually 

<*V> « (^V)iatt - ^e^, (8) 
where we have subtracted the dominant ultraviolet divergence. 
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With these approximations, the dynamics of the system depends only 
on two unknown parameters, the Higgs mass to#, for which we used the 
value rriH = 100 GeV, and the initial value fa, which is constrained by the 
requirement that when the system equilibrates, it is already deep enough in the 
broken phase to prevent the washout of the baryon asymmetry by sphalerons. 
It has been estimatccB that this requires 0>T ro h C atj where T ro h at is the final 
temperature. On the other hand, conservation of energy implies 



Since we were only interested in the qualitative behaviour and not in precise 
numbers, we used in the simulations the value fa = 700 GeV, which is slightly 
larger than the constraint (|l(]) . 

3 Results 

Startng from the initial configuration described above, we solved numerically 
the equations of motion (|l|) on a 60 3 lattice with lattice spacing Sx = 3 TeV -1 
and time step St = 0.6 TeV -1 . The time evolution of \4>\ 2 is shown in Fig. |[ Its 
qualitative behaviour is similar to that in the Abelian theory.B Note that in the 
presence of fluctuations, |0| 2 is never zero, but from its qualitative behaviour 
we can deduce that the clcctrowcak symmetry is effectively restored until t ~ 
0.8 GeV -1 . In the absence of a Higgs condensate, the damping term is namely 
expected to cause |</>| 2 to decrease as exp(— Ft), which is exactly what we 
observe. At t ~ 0.7 GeV -1 , \(f>\ 2 starts to grow towards its vacuum expectation 
value, just as it is expected to do in the broken phase. During this period of 
non-thermal symmetry restoration, baryon number is not conserved, and the 
out-of-equilibrium processes can generate a non-zero baryon asymmetry. 

4 Baryon asymmetry 

So far, we have concentrated on understanding the qualitative dynamics of 
the clcctrowcak theory during preheating. However, IfjKC really want to test 
the scenario of electroweak baryogenesis at preheatingj3'Q we have to be able 
to measure the baryon asymmetry generated during the transition, and that 
involves many technical problems. 
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which leads to the constraint 



fa < 600 GeV. 



(10) 
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Figure 1. The time evolution of \<f>\ 2 with the initial value ij>0 = 700 GeV. The dashed 
line shows the vacuum expectation value. Until t ~ 0.8 GeV -1 , the curve decreases ex- 
ponentially, indicating that the Higgs condensate is absent and the symmetry is restored. 
Eventually, |0| 2 starts to grow towards it vacuum value, which means that the condensate 
develops and the symmetry is broken. 



In principle, we can measure the change of the baryon number even though 
we don't have fcrmions in our system, because a quantum anomaly links it to 
the changes of the Chern-Simons number of the SU(2) gauge field 

AB = 3AN CS = -L^J^dtJ d 3 xe ijk E?F°- k . (11) 

By measuring Ef and .F? fc , we could then find the change of the baryon num- 
ber. 

However, as -/Vcs is a topologi caL quantity, it does not have a natural defi- 
nition on a lattice, and attemptsEJO to measure tijkEfF^ in lattice theories 
have shown that it is dominated by ultraviolet fluctuations. However, at least 
in thermal equilibrium, it is possible to remove these fluctuations by cooling 
the system,EJ which leads to a more reliable result. 

Another, more serious problem is that in order to generate the observed 
baryon asymmetry, the theory must violate CP. If this effect arises from heavy 
degrees of freedom, it can be approximated by an effective term in the La- 
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grangian 



<5cp 




(12) 



where M is the mass of the heavy fields and Sep parameterizes the strength 
of the CP violation. Unfortunately, it is very difficult to add this term in the 
equations of motion, because it is extremely sensitive to ultraviolet fluctua- 
tions. 

Thus, it seems that the only way to measure the generated baryon asym- 
metry is to treat Sep as a linear perturbation. One way to do that is to use 



and measure r sp h, T e g and (0 2 ) in the simulation. However, this approxima- 
tion may not always take all the relevant effects into account S3 

5 Conclusions 

We have studied numerically some aspects of the behaviour of the electroweak 
theory during preheating and found that it is possible to restore the symme- 
try non-thermally for a short time, which allows the baryon asymmetry to 
be generated. When the Higgs and gauge bosons decay into fermions, the 
temperature decreases so rapidly that this baryon asymmetry does not have 
time to be washed out. 

While the results presented in this talk support the scenario of electroweak 
baryogencsis at preheating, they arc qualitative in nature and do not let us 
deduce the generated amoixnt of baryon asymmetry. However, more precise 
simulations are under wayO 
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